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Abstract 

This paper describes a study of a lean premixed (LP) methane-air combustion 
wave in a two-dimensional Cartesian and axisymmetric coordinate system. Lean 
premixed combustors provide low emission and high efficiency; however, they are 
susceptible to combustion instabilities. The present study focuses on the behavior of the 
flame as it interacts with an external acoustic disturbance. It was found that the flame 
oscillations increase as the disturbance amplitude is increased. Furthermore, when the 
frequency of the disturbance is at resonance with a chamber frequency, the instabilities 
increase. For the axisymmetric geometry, the flame is found to be more unstable 
compared to the Cartesian case. In some cases, these instabilities were severe and led to 
flame extinction. In the axisymmetric case, several passive control devices were tested to 
assess their effectiveness. It is found that an acoustic cavity is better able at controlling 
the pressure fluctuations in the chamber. 


Introduction 

Environmental awareness leads to increased demands for the ultra-low emission 
engines. Lean premixed (LP) combustion is the most effective process for lowering the 
emission. As a result, the research and development of gas turbines, which use LP 
combustors, have been a major interest for the past two decades. Future gas turbine 
engines must provide ultra-low emissions and high efficiency at low cost while 
maintaining the reliability and operability of present day engines. The demands for 
increased performance and decreased emissions have resulted in advanced combustor 
designs that are critically dependent on efficient fuel/air mixing and lean operation. 
However, all combustors, but most notably lean-burning low-emissions combustors, are 
susceptible to combustion instabilities. These instabilities typically occur due to sustained 
pressure perturbations of acoustic nature in a chamber where combustion takes place. 
They arise when different chamber processes couple in such a way that inherent 
disturbances in the system are self-excited via interactions with the combustion process. 
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The acoustic pressure alters the flame, which in turn drives sound by producing unsteady 
heat release. 


Frendi [1] concentrated on the coupling between the unsteady heat release rate 
from combustion and the acoustic energy of the combustion chamber. Furthermore, he 
also studied the effect of adding a Helmholtz resonator to the system. Similarly, Paxson 
and Quinn [2] performed a comprehensive study of the effect of pressure waves on the 
acoustic modes of the combustion chamber and how it instigates various forms of 
instabilities. In addition, Paxson [3] studied the effect of abrupt change of area in a lean 
premixed combustor rig and how this change of area expedites the phenomena of 
combustion instabilities. On the other hand, Prasad [4] investigated the effect of external 
pressure perturbations on “premixed flames” using detailed chemistry. His results 
indicated that reflected and transmitted waves differed significantly from incident waves. 
However, Wu, et. al. [5] studied the flame structure and its associated hydrodynamic field 
in detail. They proposed that combustion instability arises when acoustic modes of the 
chamber are excited and amplified by the flame instabilities through mutual resonance. 

In the present study, a two-dimensional Cartesian and axisymmetric model is 
used to study the effects acoustic coupling on combustion instabilities. A one-step global 
reaction is used to model the chemical kinetics. 


Mathematical Model 

The Mathematical model consist of the nonlinear Navier-Stokes equations with a 
one-step irreversible exothermic chemical reaction of the form 


v F (Fuel) +v 0 (Air) ->v p (Products) 


( 1 ) 


where Fuel and Air are assumed to be in gas phase. 


A set of unsteady, compressible Navier-Stokes equations for a two-dimensional 
axisymmetric system used to describe our gas-phase reacting flow problem, is as follows 


23.+*+*. B 

dt dx dr 


( 2 ) 


where Q is the state vector, F and G are the fluxes in the x and r directions, respectively, 
and B is the source term that arises in cylindrical coordinates in addition to the reacting 
flow. 




pu 


pv 



’ p " 

pu 

pv 

II 

P“ 2 -*xx+P 

puv~r xr 

u(pe+p)-l— 

G = r 

puv-T xr 

P» 2 -*rr+P 

v(pe+p)-A^- 

B = 

0 

0 

P~ x ee 

pe 

PY F 


ox 

v ^y f 

puYp-pD-rr 



rqcoJ”+rq ig " 

4 Ml 

rniF 

pY 0 




pvY 0 -pD~f- 

or 


* Ml 

rmo 



( 3 ) 


where, Y 0 and Y F are Oxygen and fuel mass fractions, respectively. The source terms in 
the energy equations are given by 
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Results and Discussions 

The System of equations given by (2) - (3) with source terms (4) - (7) is solved 
using the Gottlieb-Turkel scheme [6], which is second order accurate in time and fourth 
order accurate in space. This scheme is based on MacCormack second order accurate 
Scheme. Along with Gottlieb-Turkel scheme, fourth order artificial dissipation terms 
given by Anthony Jameson [7] were also explicitly added to capture sudden 
discontinuities. To solve the above system of equations, coding is done in Matlab. A 
domain having of 6 cm in each direction is used with 200x200 grid points. At first the 
acoustic modes of the 2D Cartesian as well as axisymmetric chamber are determined 
using a numerical bomb test. This is achieved by using the source term given by equation 
(5). In this numerical experiment, the combustion source terms are ignored. The modes of 
the chamber are then obtained by taking the FFT of the time history of the pressure time 
history. To validate our work, a special case was studied using an axisymmetric chamber 
with a diameter and length of 2.125 and 8.5 inches, respectively. The new chamber 
geometry corresponds to the chamber used by Moser et al. [8] in their experiments. 

Validation Case: 


Figure 1 shows the temperature contours in the axisymmetric chamber shortly after 
ignition. A symmetrical high temperature kernel is shown at the center of the chamber. 
Figure 2 shows the corresponding instantaneous pressure contours. The contours show 
clearly the presence of a longitudinal mode. This is confirmed by the time history and the 
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corresponding power spectral density of Figure (3a&b) showing a strong mode at around 
3000 Hz. These results agree well with those obtained experimentally by Moser et al. [8], 
Fig (3 c), though our frequency is slightly higher because of a small difference in the 
chamber dimensions. 

Combustion-Acoustic Interaction Results: 


Figures (4)a-b and (5)a-b show the pressure time history at the center of the chamber and 
the corresponding FFT for both the Cartesian (Fig.4(a)-(b)) and axisymmetric (Fig. 5(a)- 
(b)) geometries, respectively. The pressure signature shown in figures (4) and (5) is 
obtained from a numerical bomb experiment. The time histories show the presence of 
several frequencies which is confirmed by the power spectral densities. The 
corresponding temperature history at the center of the Cartesian chambers are shown on 
Fig. (6). The figure shows a transient behavior of the temperature. The instantaneous 
temperature contours corresponding to the Cartesian geometry is shown on Fig. (7) for 
the case without acoustic disturbance. From Fig.(7), it is observed that a symmetric 
kernel is obtained. When an acoustic disturbance is interacting with the combustion 
wave, depending upon the amplitude and the frequency of the disturbance, it is found that 
the flame is strongly affected by the acoustic disturbance as shown on Fig. (8) and (9). 


Passive Control 


Figure 10 shows a composite plot of the grid used, the pressure contours, the time history 
of the pressure and die corresponding power spectral density for an axisymmetric 
combustion chamber without a passive device. A strong fundamental and harmonic are 
shown on the power spectral density. When a quarter-wave is added to the chamber, 
Figure 1 1 shows that the pressure fluctuations are significantly reduced and the power 
spectrum shows a large number of peaks with a weakened fundamental and harmonic. 
Figure 12 shows the effect of using an acoustic cavity on the resonant interaction. Further 
reduction in the amplitude of the pressure fluctuations is obtained. Both passive control 
devices act to detune the chamber. 


Conclusions 

The results show that acoustic-combustion coupling plays an important role in 
combustion instability. The response of a premixed flame as it interacts with an 
externally imposed acoustic disturbance were studied. Based on the results presented, the 
following conclusions are made: 

1 . The instabilities of the flame increase by increasing the disturbance amplitude. 

2. When the frequency of the disturbance is at resonance with the chamber 
frequency, the instabilities increase and lead to flame extinction. 

3. For the axisymmetric geometry, the flame behaves more unstably compared to the 
one for the Cartesian system. 
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4. An acoustic cavity is a more effective passive control device than either a baffle 
or quarter wave. 

Comparisons to experimental results show reasonable agreement. 
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Figure 1: Instantaneous temperature contours of a premixed-flame in a 2-D axisymmetric chamber 

at t = 0.02 seconds. 
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Figure 2: Instantaneous pressure contours of a premixed-flame in a 2-D axisymmetric chamber 

at t = 0.02 seconds. 
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Figure 3: (a) Time history of the pressure of a lean-premixed flame in a 2-D 
axisymmetric geometry, (b) the FFT of the corresponding pressure data for the special 
case and (c) The frequency vs. time plot from the experiment results by Moser [8] 




(b) 


Figure 4: (a) Time history of the pressure of a lean-premixed flame in a two-dimensional 
Cartesian geometry and (b) the FFT of the corresponding pressure data 




Figure 5: (a) Time history of the pressure of a lean-premixed flame in a two-dimensional 
axisymmetric geometry and (b) the FFT of the corresponding pressure data 
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Figure 6: Time history of the temperature of a lean-premixed flame in a two-dimensional 
Cartesian geometry 



Figure 7: Instantaneous temperature contours of a premixed-flame with no acoustic disturbance at 
t = 0.005 seconds for the Cartesian geometry 
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(C) 

Figure 8: Instantaneous temperature contours of a combustion-acoustic interaction in a 2-D 
Cartesian system at t = 0.005 seconds for (a) s= 0 , (b) e=0.3 and (c) s=0.6. 



(c) 


Figure 9 - Instantaneous temperature contours of a combustion-acoustic interaction in a 2-D 
Cartesian system at t = 0.0025 seconds for (a) no-disturbance , (b) off-resonant and (c) resonant case 
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Figure 10: Resonant acoustic wave interaction with a combustion chamber mode. 
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Figure 1 1 : Effect of a quarter- wave on the resonant interaction between an acoustic wave 
and a chamber mode. 
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Figure 12: Effect of an acoustic cavity on the resonant interaction between an acoustic 
wave and a chamber mode. 
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